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The channeling effect in crystals refers to the orientation dependence of charged ion penetration in crystals. 
In direct dark matter crystalline detectors, a channeled ion recoiling after a collision with a WIMP gives all 
its energy to electrons. Thus channeling increases the ionization or scintillation signal expected from a WIMP. 
Channeling is a directional effect which depends on the velocity distribution of WIMPs in the dark halo of our 
Galaxy and could lead to a daily modulation of the signal. I will present estimates of the expected amplitude 
of the daily modulation in direct dark matter detectors, both due to channeling and just due to the rotational 
velocity of the Earth around itself. 

o 

O-t, 1. Introduction 

D 

Understanding ion channeling and blocking in crystals is important in the interpretation of direct dark matter 
experiments measuring ionization or scintillation signals. In the "channeling effect" ions incident upon a crystal 
along symmetry axes and planes suffer a series of small-angle scatterings that maintain them in the open 
"channels" between the rows or planes of lattice atoms and thus lose their energy predominantly to electrons 
(while non-channeled ions transfer most of their energy to lattice nuclei). The "blocking effect" consists in a 
reduction of the flux of ions originating in lattice sites along symmetry axes and planes due to the shadowing 
effect of the lattice atoms directly in front of the emitting lattice site [l| . 

The potential importance of the channeling effect for direct dark matter detection was first pointed out for 
Nal (Tl) by Drobyshevski Q and subsequently by the DAMA collaboration Q. When Na or I ions recoiling 
<— » after a collision with a dark matter WIMP (Weakly Interacting Massive Particle) move along crystal axes and 
planes, their quenching factor is approximately Q = 1 instead of QNa = 0.3 and Q\ — 0.09, since they give 
their energy to electrons. The DAMA collaboration estimated the fraction of channeled recoils and found it to 
be large for low recoiling energies in the keV range. This effect considerably shifts the regions in cross section 
versus mass of acceptable WIMP models in agreement with the DAMA data towards lower WIMP masses. As 
a consequence, the comparison between the DAMA results and the null results of other experiments is affected 
too 01 . 

• In 2008, Avignone, Creswick, and Nussinov Q suggested that ion channeling in Nal crystals could give rise 
to a daily modulation of the dark matter signal. Such a modulation of the rate due to channeling is expected 
. to occur at some level because the "WIMP wind" arrives to Earth on average from a particular direction fixed 
to the Galaxy. Assuming that the dark matter halo is on average at rest with respect to the Galaxy, this is the 
direction towards which the Earth moves with respect to the Galaxy. Earth's daily rotation naturally changes 
the direction of the "WIMP wind" with respect to the crystal axes, thus changing the amount of recoiling ions 
that are channeled vs. non-channeled. This amounts to a daily modulation of the dark matter signal detectable 
■ via scintillation or ionization. In addition to the daily modulation due to channeling, there exists a daily 
modulation in the usual (i.e. non-channeled) signal rate which is due purely to the change of WIMP kinetic 
energy in the lab frame as the Earth rotates around itself. This means that the number and energy of WIMPs 
above threshold in the lab frame changes during a day. If measured, the daily modulation due to channeling or 
in the usual rate would be a background free dark matter signature. 

My collaborators, Graciela Gelmini and Paolo Gondolo, and I @ have evaluated the upper bounds on the 
channeling fractions for different crystalline detectors used in dark matter experiments through analytical means. 
To calculate the amplitude of the daily modulation in Nal crystals, we compute the event rate during a day 
using the channeling fractions and the actual differential recoil spectrum and taking into account channeled and 
non-channeled recoils. We also examine the possibility that such a daily modulation might be observable in the 
data accumulated by the DAMA collaboration. 



> 

in 
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2. Channeling Fractions 

Our calculation of the channeling fractions as function of recoil energy is based on the classical analytic 
models which started to be developed in the 1960's and 70's. In particular we use Lindhard's model Q, 
supplemented by the planar channel models of Morgan and Van Vliet H and the work of Hobler Q on low 
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Figure 1: (a) Channeling fraction for a 50 keV Na ion in different directions plotted on a sphere using HEALPix: 
probability equal to one in red, and probability equal to zero in blue, (b) Fraction of channeled incident I (black) and 
Na (green/gray) ions as a function of their incident energy E with the static lattice without (dot dashed lines) and with 
(solid lines) dechanneling due to interactions with Tl impurities. The results of DAMA are also included (dashed lines). 



energy channeling. In these models the rows and planes of lattice atoms are replaced by continuum strings and 
planes, in which the screened Thomas-Fermi potential is averaged over a direction parallel to a row or a plane. 
The continuous potential, U is considered to be uniformly smeared along the row or plane of atoms, which 
is a good approximation if the propagating ion interacts by a correlated series of many consecutive glancing 
collisions with lattice atoms in the row or plane. Just one row or plane is considered in this model. Lindhard 
proved that for an ion propagating with kinetic energy E, and for small angle <j> between the ion's trajectory 
and the atomic row (or plane) in the direction perpendicular to the row (or plane), the so called "transverse 
energy" , E±_ = E sin 2 (j> + U ~ E(\? + U is conserved. 

The conservation of the transverse energy provides a definition of the minimum distance of approach to the 
string (or plane) of atoms, p m i n , at which the trajectory of the ion makes a zero angle with the string (or plane), 
and also of the angle ip at which the ion exits from the string (or plane), i.e. far away from it where U ~ 0. 
Channeling requires that p m m > Pc(E), where p c (E) is the smallest possible minimum distance of approach of 
the propagating ion with the row (or plane) for a given energy E. This amounts to ip < ip c (E), where ipc(E) is 
the maximum angle the ion can make with the string far away from it (i.e. in the middle of the channel) if the 
ion is channeled. 

2.1. Channeling of Incoming Particles 

The channeling of ions in a crystal depends not only on the angle their initial trajectory makes with strings 
or planes in the crystal, but also on their initial position. Ions which start their motion close to the center of 
a channel, far from a string or plane, where they make an angle ip, are channeled if the angle is smaller than 
a critical angle (as explain earlier) and are not channeled otherwise. Particles which start their motion in the 
middle of a channel (as opposed to a lattice site) must be incident upon the crystal. 

No data or simulations of Na and I ions propagating in a Nal crystal is available at low energies. We show 
that to a good approximation we can use analytic calculations and reproduce the channeling fraction in Nal 
presented in Ref. [3| . Temperature corrections are neglected and we use a static lattice (similar to the approach 
of Ref. Q). For an incident angle ip with respect to each of the channels and an ion energy E, the fraction 
of channeled incident ions for axial and planar channels is 1 if ip is smaller than the critical angle for the 
corresponding channel and zero otherwise. 

To find the total fraction of channeled incident nuclei, we integrate the direction-dependent channeling fraction 
over the incident direction. The integral cannot be solved analytically, so we integrate numerically by performing 
a Riemann sum once the sphere of directions has been divided using the Hierarchical Equal Area iso-Latitude 
Pixelization (HEALPix) [lCj ]. HEALPix uses an algorithm to deal with the pixelization of data on a sphere, and 
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it is useful to compute integrals over direction by dividing the surface of a sphere into many pixels, computing 
the integrand at each pixel (i.e. each direction, see Fig. [TJa), and finally summing up the values of all pixels 
over the sphere. 

A channeled ion can be pushed out of a channel by an interaction with an impurity such as the atoms of Tl 
in Nal (Tl). Here we simply assume that if a channeled ion interacts with a Tl atom it becomes dechanneled 
and thus it does not contribute to the fully channeled fraction any longer. We thus neglect the possibility that 
after the interaction the ion may reenter into a channel, either the same or another. Figure [TJa shows the axial 
and planar channels of the Nal crystal for incoming Na ions with an energy of 50 keV. We include here only the 
channels with lower crystallographic indices, i.e. 100, 110 and 111, which provide the dominant contribution 
to the channeling fraction, as is also done in Ref. Q. The fraction of channeled incident Na or I ions as a 
function of their incident energy E is shown in Fig. [TJb. For comparison, Fig. [TJb also shows the channeling 
fraction obtained by DAMA. Good agreement with the channeling fractions of DAMA is achieved only when 
dcchanneling due to the interaction with Tl impurities is included. 



2.2. Channeling of Recoiling Lattice Nuclei 

The recoiling nuclei start initially from lattice sites (or very close to them), thus blocking effects are important. 
In fact, as argued originally by Lindhard [7], in a perfect lattice and in the absence of energy- loss processes the 
probability that a particle starting from a lattice site is channeled would be zero. The argument uses statistical 
mechanics in which the probability of particle paths related by time-reversal is the same. Thus the probability 
of an incoming ion to have a particular path within the crystal is the same as the probability of the same ion to 
move backwards along the same path [l| . This is what Lindhard called the "Rule of Reversibility" . Using this 
rule, since the probability of an incoming channeled ion to get very close to a lattice site is zero, the probability 
of the same ion to move in the time-reversed path, starting at a nuclear site and ending inside a channel is also 
zero. However, any departure of the actual lattice from a perfect lattice, for example due to vibrations of the 
atoms in the lattice, would violate the conditions of this argument and allow for some of the recoiling lattice 
nuclei to be channeled, as already understood in the 70's (ll| . 




In our model, a recoiling ion is channeled if the collision ion- WIMP happens at a distance large enough from 
the string or plane to which the ion belongs. Namely, channeling happens if the initial position of the recoiling 
motion is pi > pi, min - We define E± in terms of the initial recoil energy E of the propagating ion, the angle of 
the initial recoil momentum with respect to the particular string or plane of atoms 4>i, and the initial position 
Pi. Since the potential U{p) decreases monotonically with increasing p, U(p m i n ) < U(p c (E)), and the condition 



i.e. E±(E,4>i, pi, min) = U{p c {E)) defines the distance Pi jm i n in terms of p c . 

We take the initial distance distribution, pi of the colliding atom to be a Gaussian with a one dimensional rms 
vibration amplitude u\ (using the Debye model), and we obtain the probability of channeling for each individual 
channel by integrating the Gaussian between the minimum initial distance and infinity (a good approximation to 
the radius of the channel) . The dependence of these probabilities on the critical distances enter in the argument 
of an exponential (for axial channels) or an erfc function (for planar channels). In order to obtain the total 
geometric channeling fraction we sum over all the individual channels we consider. The integral over initial 
directions is computed using HEALPix. 

Figure [2] shows what we consider to be our main predictions for the range expected as an upper limit to 
the channeling fraction in Nal at 239 K. The parameter c mentioned in the figure is a number that we expect 
to be between 1 and 2, which regulates the importance of temperature corrections. The channeling fractions 
are typically smaller for larger values of c. Notice that we have not included any dechanneling effects due to 
impurities or dopants, which should also decrease the channeling fractions. Using the fractions in Fig. [5J the 
difference in the regions of cross section versus mass of acceptable WIMP models corresponding to the DAMA 
annual modulation signal differs only at la if channeling is included or not included In the absence of 
channeling, the DAMA region could be compatible with the CoGeNT region for light WIMPs 12]. 

Our results for the total geometric channeling fraction for Si ions propagating in a Si crystal and Ge ions 
propagating in a Ge crystal at different temperatures are shown in Figs. [3] and 0] for the two cases of c = 1 and 
c = 2, respectively. The channeling fractions we obtain are strongly temperature dependent. As the temperature 
increases, the probability of finding atoms far from their equilibrium lattice sites increases, which increases the 
channeling fractions, but the critical distances p c become larger which decreases the channeling fractions. 




E ± (E,(f>i,pi) = U(p min ) < U(p c ) 



(1) 
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Figure 2: Upper bounds to the channeling fraction of recoiling Na (solid lines) and I (dashed lines) ions in Nal for c = 1 
(black) and c = 2 (green/gray) cases at T=293 K. 
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Figure 3: Channeling fractions of (a) Si and (b) Ge recoils in a Si and a Ge crystal respectively, as a function of the ion 
energy for temperatures T = 900 °C (orange or medium gray), 600 °C (green or light gray), 293 K (black), and 44 mK 
(blue or dark gray) in the approximation of c = 1. 
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Figure 4: Same as Fig. [3] but for c = 2. 
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3. Angular Distribution of Recoil Directions due to WIMPs 



Consider the WIMP-nucleus elastic collision for a WIMP of mass m and a nucleus of mass M. The 3- 



spectrum as a function of the recoil momentum 



dR 



elocity 

q[H 



dE R dQ q 



ATTlTlfl 2 a \2/i' 



(2) 



where E R is the recoil energy, dfl q = d(f>d cos 9 denotes an infinitesimal solid angle around the recoil direction 
q = q/q, q = |q| is the magnitude of the recoil momentum, [i = mM/(m + M) is the reduced WIMP-nucleus 
mass, q/2fi is the minimum velocity a WIMP must have to impart a recoil momentum q to the nucleus, p is 
the dark matter density in the solar neighborhood, do is the total scattering cross section of the WIMP with 
a (fictitious) point-like nucleus, and S(q) is the nuclear form factor normalized to 1. We concentrate here on 
WIMPs with spin-independent interactions. 

For a truncated Maxwellian WIMP velocity distribution with respect to the Galaxy, with dispersion a v and 
truncated at the escape speed v csc , the Radon-transform is [l3| 



/lab[ ^,qj - 



1 



iY eS c(27ra2)i/2 



exp 



[(q/2n) + q.V 
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if (q/2/j,) + q.Vi a b < v csc , and zero otherwise, where 



N e , a = erf 



\V2~<Jv 



2 v e 



7r a v 



■ exp 



^csc 
2^2 



(4) 



Here we are assuming the detector has a velocity Vi a b with respect to the Galaxy (thus — Vi a b is the average 
velocity of the WIMPs with respect to the detector). Vi a b is defined in terms of the galactic rotation velocity 
VcaiRot at the position of the Sun (or Local Standard of Rest (LSR) velocity) , Sun's peculiar velocity Vsoiar in 
the LSR, Earth's translational velocity VsarthRev with respect to the Sun, and the velocity of Earth's rotation 
around itself V Earth Rot- Thus, Vi ab = V Ga iRot + V So iar + V Ear thRov + VF.art.hRnt.. We take V Ga iRot either 220 
km/s or 280 km/s, as reasonable low and high values (as done in Ref [14]]), which correspond to T4 a b either 
228.4 km/s or 288.3 km/s, respectively. The presence of q.V\ a b in Eq. [3] means that in order to compute the 
differential rate we need to orient the nuclear recoil direction q with respect to V"i a b- 
The differential energy spectrum is given by 



dR 
dE 



dR 



dEjidflq 



P {E,E R ,q)dn q dE R , 



(5) 



where p(E, E R , q)dE is the probability that an energy E is measured when a nucleus recoils in the direction q 
with initial energy Er. With our analytic approach we cannot estimate the importance of dechanneling mech- 
anisms, such as the presence of lattice imperfections, impurities or dopants. Thus we disregard dechanneling, 
and assume that a recoiling nucleus can only either be channeled, in which case the measured energy is the 
whole initial recoil energy E = E R (first term in the following equation) or not channeled, in which case the 
measured energy is E = QEr (second term), 



p(E, E R , q) = X (E R , q)6(E - E R ) + [1 - X (Er, q)]<5(£ - QE R ), 



(6) 



where x(E R , q) is the probability that a nucleus with recoil energy E R is channeled in a given direction q. 

Inserting Eqs. [2] and [5] in Eq. [SJ we find the measured differential rate. The integral in Eq. [5] cannot be 
computed analytically, and we integrate numerically using HEALPix. 



4. Daily Modulation in Nal Crystals 

We present here the daily modulation amplitude expected in Nal crystals, assu ming that WIMPs have a 
truncated Maxwellian velocity distribution with i> csc = 650 km/s and o~ v = 300 km/s [13J. 
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Figure 5: Signal rate (in events/kg-day-keVee) as function of the Universal Time (UT) during 24 hours for m — 10 GeV 
and (a) E = 10 keVee and Viab = 228.4 km/s, and (b) E = 6 keVee and Viab = 288.3 km/s. The parameters used are 



= 0.2, Qj = 0.09, <j„ = 2 x 10" 



c = 1 for temperature effects and a crystal temperature of T = 293 K. 



The spin- independent detection rate of WIMPs has a time dependence through the Radon transform /i a b. 
Notice that /i a b (see Eq. [3]) changes during a day through the (q.Vi a b) factor appearing in the exponent and 
the dependence of Vi ab on V Ear thRot- 

Here we show the signal rate as a function of time during a particular arbitrary Solar day (September 25, 
2010). The total rate consists of signal plus background, Rt = R s + Rb, and we assume that there is no daily 
modulation in the background. We define the relative signal modulation amplitude A s (taking into account the 
signal only) in terms of the maximum and minimum daily signal rate R s as 

* ^s-max Rs— min /~\ 

S= R 7r ' ^ ' 

rt s-max ~r -fi-s — min 

The average rate due to the signal alone is R s — (R s -max + #s-min)/2. 

Exploring the parameter space of WIMP mass and WIMP-proton cross section for different recoil energies 
we find that the relative modulation amplitudes A s can be large, even more than 10% for some combination 
of parameters. We show two examples in Fig. [5J where we plot the signal rate (in events/kg/day/keVee) as a 
function of the Universal Time (UT) during 24 hours. We find that the largest A s happen when the signal is 
only due to channeling. This happens when there are no WIMPs in the galactic halo with large enough kinetic 
energy to provide the observed energy if the recoil is not channeled. 

The detectability of a particular amplitude of daily modulation depends on the exposure and background 
of a particular experiment. The former DAMA/Nal and the DAMA/LIBRA experiments have a very large 
cumulative exposure, 1.17 ton x year. However even with this large exposure, we find that the daily modulations 
we predict are not observable. To observe the daily modulation, the total number of events Nt (N s signal plus 
iVj, background events) over the duration of the experiment should be divided into two bins, the "high-rate" 
bin with Nx-max events and the "low-rate" bin with Nx-min events, so that Nt — Nt— max + iVr-min' F° r the 
daily modulation to be observable at the 3u level one should have 



A T _ max " N T -mm = A S N S > 3(7 ~ 3y/N T /2, (8) 

where a 2 ~ A^/2 because, with a small modulation, on average A^-max — iVr-min — At/2. 

If the detector exposure is MT in kg-day and we take bins of width AE in keVee, then At = RtMT AE 
and N s = R S MTAE, where the rates are in events/kg-day-keVee. Thus the condition in Eq. [5] becomes 

R 2 JR T > 9/(2A 2 s MT AE). (9) 

The total rate of the DAMA experiment at low energies between 4 keVee and 10 keVee is Rt — 
1 events/kg/day/keVee [HI]. This rate is much larger than the signal rates we predict and is, therefore, dom- 
inated by background. We choose here an energy bin AE ~ 1 keVee, narrow enough to assume the signal 
rate to be constant in it and compatible with the energy resolution of DAMA [l6j . With the cumulative 
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exposure of DAMA, the condition in Eq. [§] for relative daily modulation amplitude A s observable at 3<r is 
R s A s > 3.2 x 1CT 3 events/kg/day/keVee or 

-R s -max - Rs-min > 6.4 x 1CT 3 events/kg/day/keVee. (10) 

Even the largest relative daily modulations we find, shown in Fig. [5J are not observable in the DAMA data 
according to Eq. [TU] 

One could ask what exposure would be needed with the current total rate in the DAMA experiment to make 
the daily modulation observable. Setting Rt — 1 events/kg/day/keVee in Eq. [9l we obtain 

MTAE 9 18 

(events/kg/day/keVee) > 2 (A s R s ) 2 ~ (fl s _ max - i? s - min ) 2 ' 

For the case with the highest rate difference we found (to = 10 GeV, E = 6 keVee and Vj a b = 288.3 km/s, shown 
in Fig. [SJb) and with Ai? ~ 1 keVee we would require an exposure 40 times larger than the current exposure of 
DAMA. 

5. Future Prospects for Other Experiments 

We have computed the daily modulation due to channeling in other material such as Ge, solid Xe and solid 
Ne, and we find that it will be very difficult to observe. For light WIMPs the cross section can be larger than 
for heavier ones without violating experimental bounds, a p = 10 _39 cm 2 [13, EH and this favors the detection of 
the daily modulation. We find that for a WIMP mass m — 5 GeV the daily modulation due to channeling may 
be observable in solid Ne if the signal would be above threshold and assuming no background. For example for 
a solid Ne detector operating at 23 K at Gran Sasso, for E = 10 keV, assuming Qnc = 0.25 [T^|, c = 1 and 
Vi a b = 228.4 km/s, we find that the exposure needed to observe the modulation at 3cr is MT = 0.33 ton year. 
The usual rate is zero in this case, and the modulation is just due to channeling. 

The daily modulation in the usual rate may be possible to observe in future non-directional detectors with 
negligible background. The required exposures to have an observable daily modulation at the 3cr level in future 
solid or liquid detectors are shown in Table Q] for cases in which channeling is negligible (or absent in liquid 
detectors), and the daily modulation is due purely to the change of WIMP momentum distribution in the 
lab frame as the Earth rotates around itself. The exposures are calculated assuming no background and with 
Vi a b = 228.4 km/s for a Ge detector at the Soudan Underground Laboratory (with Qq c — 0.2 |17fl) and all other 
detectors at Gran Sasso. 



Table I: Observability in Future Detectors 



Detector 


E R (keV) 


m (GeV) 


a v (cm 2 ) 


MT (ton year) 


Ge 


10 


7 


1Q -41 


11 


Solid Ne 


8 


5 


1Q -39 


0.8 


Liquid Ne 


14 


5 


1Q -39 


3.4 


Liquid Xe 


3 


5 


1Q -39 


0.09 


Liquid Ar 


8 


5 


1Q -39 


0.9 



6. Conclusions 

We have studied the channeling of ions recoiling after collisions with WIMPs in different crystalline detectors 
such as Nal (Tl), Si and Ge. Channeled ions move within the crystal along symmetry axes and planes and 
suffer a series of small-angle scatterings that maintain them in the open "channels" between the rows or planes 
of lattice atoms and thus penetrate much further into the crystal than in other directions. Ions which start 
their motion close to the center of a channel, at an initial angle ip, are channeled if the initial angle is smaller 
than a critical angle, and are not channeled otherwise. With the simple model of dechanneling we used for 
Nal, we could reproduce the channeling fractions computed by the DAMA collaboration. However this is not 
the case in direct dark matter detectors, since the recoiling ions start their motion at or close to their original 
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lattice sites. For recoiling ions, blocking effects (neglected in the DAMA calculation) are important, and the 
channeling fraction is smaller. As argued originally by Lindhard, in a perfect lattice and in the absence of 
energy-loss processes, the probability that a particle starting from a lattice site is channeled would be zero. 
However, due to vibrations in the crystal, the atom that interacts with a WIMP may be displaced from its 
position in a perfect lattice, and there is a non-zero probability of channeling. As seen in Fig. [5] for Nal, without 
including dechanneling, the channeling fraction is never larger than 5%. 

We have also studied the possibility of a daily modulation due to channeling or in the usual non-channeled 
rate, which could be a background free signature of dark matter. We find large daily modulation amplitudes 
for Nal (even more than 10%) which are not observable in DAMA. However, the daily modulation might be 
detectable in other experiments such as Ge, solid and liquid Ne, and liquid Xe and Ar with smaller background 
or larger exposure. We intend to further explore the observability of a daily modulation in future experiments 
for different halo models in future work. 
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